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Abstract. We investigate the behaviour of dark energy using the recently released 
supernova data of Riess et al., 2004 and a model independent parameterization for 
dark energy (DE). We find that, if no priors are imposed on 51o m and h, DE which 
evolves with time provides a better fit to the SNe data than ACDM. This is also true 
if we include results from the WMAP CMB data. From a joint analysis of SNe+CMB, 
the best-fit DE model has wq < — 1 at the present epoch and the transition from 
deceleration to acceleration occurs at zt = 0.39±0.03. However, DE evolution becomes 
weaker if the ACDM based CMB results O 0m = 0.27 ± 0.04, h = 0.71 ± 0.06 are 
incorporated in the analysis. In this case, zt = 0.57 ±0.07. Our results also show that 
the extent of DE evolution is sensitive to the manner in which the supernova data is 
sampled. 



1. Introduction 

Supernova observations [1, 2] were the first to suggest that our universe is currently 
accelerating. Subsequently, a combination of results from cosmic microwave background 
(CMB) experiments and observations of galaxy clustering served to strengthen this world 
view [3, 4], and it is now believed that as much as 2/3 of the total density of the universe 
is in a form which has large negative pressure and which is usually referred to as dark 
energy (DE). 

The earliest theoretical model of DE - the cosmological constant (A) - satisfied 
w = p/p = — 1. Since the energy density in A does not evolve, its present value 
p A ~ 1CT 47 GeV 4 is also its initial value. As a result, the ratio Pa/ p r , where p r is the 
radiation density, had the miniscule value 10~ 123 at the Planck time. The enormous 
amount of fine tuning this might involve led theorists to suggest that, like other forms 
of matter in the universe, DE density too may show significant time-evolution. However, 
this argument for significant variability of A with redshift is questionable. Actually, it 
is a variant of the Dirac's large number hypothesis that proved not to be valid in our 
Universe. For example, the density of water has also a miniscule value of 10 -93 in Planck 
units (though, of course, larger than that of p\), and we know that its relative change 
from present time up to redshifts of the order of one (due to possible variations of the 
fine structure constant and the electron and proton masses) is less than 10~ 5 . So, based 
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on the argument above, it would be wrong to assume that dark energy density should 
significantly change with redshift. 

Another, more reliable reason to suggest a time dependent form of DE lies in the 
fact that our current accelerating epoch is unlikely to have been unique. In fact there 
is considerable evidence to suggest that the universe underwent an early inflationary 
epoch during which its expansion rapidly accelerated under the influence of an 'inflaton' 
field which, over sufficiently small time scales, had properties similar to those of a 
cosmological constant. Inspired by inflationary cosmology, quintessence models invoke 
a minimally coupled scalar field to construct a dynamically evolving model for DE. 

Recent years have seen a flurry of activity in this area and there are currently, apart 
from quintessence, at least a dozen well motivated models of an accelerating universe in 
which dark energy is a dynamically evolving quantity (see reviews [5, 6, 7]). A simple 
categorization of DE models could be as follows: 

(i) The cosmological constant, w — — 1. 

(ii) DE with w = constant ^ —1 (cosmic strings (w = —1/3), domain walls 
(w = —2/3). Quintessence with a sine hyperbolic potential [6]). 

(iii) Dynamical DE, w ^ constant. (Quintessence, Chaplygin gas [8], k-essence [9], 
braneworld models [10, 11, 12], etc.) 

(iv) DE with w < — 1. (scalar-tensor gravity models [13], phantom models, braneworld 
cosmology etc. [14, 15, 16, 17, 18].) f 

In view of the large number of possibilities for DE it would not be without advantage 
to analyse the properties of DE in a model independent manner. Such an approach was 
adopted by Alam et al.[20] (henceforth Paper I) in which the Supernova data published 
by Tonry et al.[21] and Barris et al. [22] was analysed using a versatile ansatz for the 
Hubble parameter. Paper I discovered that dynamical DE fit the SNe observations 
better than ACDM and these results found support in the subsequent analysis of other 
teams [23, 24, 25, 26]. (Paper I referred to DE evolution as 'metamorphosis' since the 
DE equation of state appeared to metamorphose from a negative present value w < — 1 
to w ~ at z ~ 1.) 

Recently Riess et al.[27] have reanalysed some earlier SNe data and also published 
new data relating to 16 type la supernovae discovered using the Hubble Space Telescope 
(HST). In the present paper we shall reconstruct the properties of DE using the new 

| Note that the kinetic energy of the phantom scalar field is negative, thus phantom models are 
strongly unstable with respect to spontaneous creation of 4 particles: a particle- antiparticle pair of this 
scalar field and a particlc-antiparticle pair of any usual matter quantum field (sec [19] for the most 
recent investigation of this process). Scalar-tensor models of dark energy with the effective Brans-Dickc 
parameter lu > —3/2 are free from this difficulty. However, a possible decrease of w below —1 in these 
models is typically of the order of 1/w, and is therefore very small, as follows from Solar System tests 
of general relativity. In this context it is interesting that the Braneworld models suggested in [11] can 
give rise to an effective equation of state for dark energy which can be substantially less than —1, and 
yet not be plagued with any instabilities at the classical or quantum level. 
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SNe data set ('Gold' in [27]). We shall also use some of the CMB results obtained by 
WMAP in the later part of our analysis. 



2. Reconstructing Dark Energy 



Perhaps the simplest route to cosmological reconstruction is through the Hubble 
parameter, which in a spatially flat universe is related to the luminosity distance quite 
simply by [28, 29, 30] 

-i -l 





We may now define the dark energy density as : 

2 

-n 0m (i + z) 3 



Pde — POc 



(1) 



(2) 



where p 0c = 3H^/(8rrG) is the present day critical density of an FRW universe, and 
fiom is the present day matter density with respect to the critical density. However, 
one should keep in mind a subtle point regarding this definition of the dark energy 
density. Its ambiguity lies in the value of Jlo m . From CMB and galaxy clustering data, 
we obtain an estimate of the total amount of clustered non-relativistic matter present 
today (denoted by fW)- However, Qo m may be different from Qo m due to a contribution 
from a part of the unclustered dark energy which also has a dust-like equation of state. 
Fortunately, this difference (if it exists) appears to be small, e.g. not exceeding 0.1 for 
the best-fit shown in section (2.2). 

Information extracted from SNe observations regarding d^z) therefore translates 
directly into knowledge of H(z), the dark energy density, and, through [31] 

d H 1 



w(x) 



aH 2 H 
2q(x) - 1 



(2x/3) H'/H 



x = 1 + z 



(4) 



3(l-On(aO) i - (H /Hyn 0m x 3 

into knowledge about the deceleration parameter of the universe and the equation of 
state of dark energy. 

For a meaningful reconstruction of DE one must construct an ansatz for H(z) which 
is sufficiently versatile to accommodate a large class of DE models. (Alternatively one 
could devise an ansatz for d^{z) or w(z); for a summary of different approaches see 
[32, 20].) An ansatz which works quite well for Quintessence and also for the Chaplygin 
gas and Braneworld models is [33] 

H(x) 



H 

= 1 — Oom 



Q 0m x d + A + Aix + A 2 x A 



X 



l + z 



(5) 



Note that this ansatz should not be considered as 



h(x) 

where A + A x + A 2 

a truncated Taylor series for h 2 (z). Rather, it is an interpolating fit for h 2 (z) having 
the right behaviour for small and large values of z. The number of terms in this fit is 
sufficient given the amount and accuracy of the present supernovae data. With more 
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and better data in future, more terms with intermediate (e.g., half-integer) powers 
of x may be added to it. Truncated Taylor series in powers of z with any finite 
number of terms may not be used to analyse most interesting quantities (including 
h(z) and di(z)) for the already existing supernovae data since these series become 
divergent at z ~ 1. In particular, for the standard ACDM model (A\ = A 2 = 0) the 
convergence radius of the Taylor expansions for h(z) and d L (z) centered at z = is 
given by the distance to the closest singular point in the complex z plane which lies 
at z = (1/fiom - l) 1/3 exp(±i7r/3) - 1. So, for Q 0m = 0.3 the Taylor series diverges at 
z > 1.197. 

This is equivalent to the following ansatz for DE density (with respect to the critical 
density) : 

Pde(x) = Pde/Poc = Aq + A x x + A 2 x 2 , (6) 

which is exact for the cosmological constant w — — 1 (Ai = A 2 = 0) and for DE models 
with w = —2/3 (A = A 2 = 0) and w = —1/3 (A — Ai — 0). Here, we neglect 
a possible difference between Q 0m and Q 0m mentioned above since it appears to be 
unimportant for the redshifts involved (z < 2), though it may become important for 
larger z. 

The corresponding expression for the equation of state of DE is : 

Aix + 2A 2 x 2 

w(x) — — 1 H — — — . (7) 

V ; 3(A + A x x + A 2 x 2 ) v ; 

A glimpse into the properties of dark energy is also provided by a two parameter 
approximation for the equation of state : 

w(z) = w + Wi z , (8) 

which can be trusted for small values of z <, 1. 

The likelihood for the parameters of the ansatz can be determined by minimising 
a x 2 -statistic: 

X (-Ho, ^Om,Pj) = 2^ 2 ' ( 9 J 

i °i 

where fio,i = itib — M = 51ogc?L + 25 is the extinction corrected distance modulus for 
SNe at redshift Zi, ai is the uncertainty in the individual distance moduli (including the 
uncertainty in galaxy redshifts due to a peculiar velocity of 400 km/s), and pj are the 
parameters of the relevant ansatz (A 1 ,A 2 for the ansatz (5) and w ,Wi for the ansatz 
(8)). We assume a flat universe for our analysis but make no further assumptions on 
the nature of dark energy. For most of the following results, we marginalise over the 
nuisance parameter H Q by integrating the probability density e~ x2 / 2 over all values of 
H . 

We first study the different subsamples of the new SNe data reported in [27] in 
some detail using the ansatz (5) and (8). Riess et al.[27] reanalysed the existing SNe 
data compiled by different search teams (mainly the High Redshift Search Team (HZT) 
and the Supernova Cosmology Project (SCP) team) and added to these 16 new SNe 




Figure 1. la, 2a, 3a confidence levels in the Wo~ w\ space for the ansatz (8) for 
fiom = 0.3, using different subsets of data from [27]. The filled circle represents 
the ACDM point. 



observed by HST. They have rejected many of the previously published SNe due to 
lack of complete photometric record, uncertain classification, etc. They divide the total 
data-set into "high-confidence" ('Gold') and "likely but not certain" ('Silver') subsets, 
fn our calculations we will consider their "high- confidence" ('Gold') subset. Figure 1 
shows the (w ,Wi) confidence levels for the ansatz (8) using different subsets of the 
new data for fl 0m = 0.3. The top two panels in fig 1 are similar to the corresponding 
panels of figure 10 of [27], and the bottom panels show other subsets of data. The panel 
(a) shows the 'Gold' sample of [27] without the 14 new HST points, and the panel (b) 
shows the full 'Gold' sample of 156 SNe. Of the last two panels, the panel (c) shows 
the results for a subset of the 'Gold' sample, consisting only of the points obtained 
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Figure 2. The logarithmic variation of dark energy density Pde(z)/po c (where 
p 0c = 3Hq/8ttG is the present critical energy density) with redshift for 
f^Om = 0.3 using different subsets of data from [27], for the ansatz (5). In each 
panel, the thick solid line shows the best-fit, the light grey contour represents 
the la confidence level, and the dark grey contour represents the 2a confidence 
level around the best-fit. The dotted line denotes matter density ^om(l + z ) 3 , 
and the dashed horizontal line denotes ACDM. 



by the SCP team and the new HST points, while the panel (d) shows the confidence 
levels for the older data set of 172 SNe published in Tonry et al. [21]. The greatest 
difference is clearly between the panels (a) and (c). We should also note that the panel 
(b) (comprising of the new reanalysed data) and panel (d) (comprising of the old HZT 
data) are qualitatively similar, with panel (b) having tighter errors. 

We now analyse these different subsets using the ansatz (5). In figures 2 and 3, we 
show the variation of dark energy density and dark energy equation of state obtained 
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Figure 3. The variation of equation of state of dark energy w(z) with redshift 
for f^om = 0-3 using different subsets of data from [27], for the ansatz (5). 
In each panel, the thick solid line shows the best-fit, the light grey contour 
represents the la confidence level, and the dark grey contour represents the 
2d confidence level around the best-fit. The dashed horizontal line denotes 
ACDM. 



using the ansatz (5) for the different subsets of data, fixing f2 0m = 0.3. It is interesting to 
note that panels (b) and (d) in both figures are quite similar, and we therefore conclude 
that the 'Gold' sample shows an evolution for DE which is consistent with that obtained 
from the older sample of Tonry et al.[21]. 

From figure 1 we find that for all four datasets, the largest degeneracy direction in 
Wq — W\ plane corresponds to the curve Wq + 0.25u>i ~ —1. This immediately suggests 
that w(z = 0.25) ~ —1, and this result appears to be quite robust since one can also 
arrive at it by choosing a very different ansatz (5) to determine w(z), as shown in 
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figure 3. 

As pointed out in [27], and confirmed by the figures 1, 2, 3, the maximum evolution 
in DE is for the 'Gold - HST' data. Less evolution is shown by the 'Gold' data set. 
However, we would like to emphasise again that the 'Gold' data set, which includes the 
14 new HST points, gives roughly the same degree of evolution for DE as the original 
data reported in Tonry et al. [21] and analysed in [20]. Thus the results pertaining to 
the evolution of DE reported in Paper I remain valid also for the new SNe (Gold) data 
set. 

Note that even the maximum rate of evolution of DE obtained in [20] in terms of 
few which is in complete agreement with the results obtained in [27]. Does 
this mean however that the new supernovae data completely exclude models with a 
fast phase transition in dark energy [34, 35] for which |g| > 1 over a narrow range 
of redshift 5z <C 1? In our view such a conclusion would be premature. The reason 
is that using an ansatz such as (5) or (8) (as in [27]) is equivalent to smoothing the 
evolution of the Universe over a redshift interval Az oc l/N, N being the number of 
free parameters used in the ansatz. Clearly, after such an implicit smoothing fast phase 
transitions in dark energy will disappear not so much because they disagree with the 
data but because of the manner in which the data has been analysed (see [20] for a more 
detailed discussion of these issues). 

Figures 1, 2, and 3 also illustrate that the degree of DE evolution can be quite 
sensitive to the manner in which the SNe data is sampled. Comparing panels (a) and 
(c) in these figures, we find that the degree of evolution of DE is largest for the 'Gold 
- HST' data set and least for the 'Gold, SCP + HST' data. (The latter is in better 
agreement with ACDM than the other three data sets.) 

2.1. Analysis of 'Gold' SNe dataset : 

We now examine the 'Gold' data set in some detail using the polynomial expansion 
of dark energy, (5). In the figure 4 we show the confidence levels for the parameters 
(Ai, A 2 ) of the ansatz for three different values of f2 0m - The x 2 value for the best-fit in 
each case is given in table 1. The \ 2 values for the corresponding ACDM models are 
also given for comparison. Interestingly, in all three cases the confidence ellipse has the 
same inclination, it only appears to shift downwards as f2 0m increases. 

In figure 5, we show the variation of the dark energy density with redshift for 
different values of the current matter density. We see that, for higher O m, the dark 
energy density evolution is sharper. The reader should also note that the growth of pde 
with time in the panels (b) and (c) is indicative of the phantom nature of DE (w < —1) 
at recent times (z < 0.25 for f2 0m = 0.3 and z < 0.4 for f2 0m = 0.4, see figure 6). 

We may obtain more information from the dark energy density by considering a 
weighted average of the equation of state: 

1 + ^ akW S [l + «»<* 111(1 + *> = IxWTT) > < 10 > 
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Figure 4. The (Al,^) parameter space for the ansatz (5) for different values 
of Oom, using the 'Gold' sample of SNe from [27]. The star in each panel 
marks the best-fit point, and the solid contours around it mark the la, 2a, 3a 
confidence levels around it. The filled circle represents the ACDM point. The 
corresponding % 2 for the best-fit points are given in table 1. 



Table 1. x 2 per degree of freedom for best-fit and ACDM models for analysis using 
the 'Gold' sample of SNe from [27]. w is the present value of the equation of state of 
dark energy in best-fit models. 

Best-fit ACDM 

^0 Xmin X 2 

0.20 -1.20 1.036 1.109 
0.30 -1.35 1.034 1.053 
0.40 -1.59 1.030 1.086 



where A denotes the total change of the variable between integration limits. Thus the 
variation in the dark energy density depicted in figure 5 is very simply related to the 
weighted average equation of state of dark energy. The value of w for different ranges 
of integration are shown in table 2. We have taken the ranges of integration to be 
approximately equally spaced in in(l + z). In all three cases shown, the value of w 
changes noticeably from close to —1 in the first bin to close to zero in the second bin. 
This indicates that the equation of state of DE is evolving from w < — 1 today to w ~ 
at z ~ 1. Note that these results are in very good agreement with those reported in 
Table 1 of Paper I. 

Figure 6 shows the corresponding variation of the DE equation of state with redshift 
for different f2 0m . Here also, there is strong evidence for evolution of DE. We see that 
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Figure 5. The logarithmic variation of dark energy density Pde(^)/Poc (where 
poc = 3Hq 1 8irG is the present critical energy density) with redshift for different 
values of f^om, using the 'Gold' sample of SNe from [27]. The reconstruction 
is done using the polynomial fit to dark energy, ansatz (5). In each panel, 
the thick solid line shows the best-fit, the light grey contour represents the la 
confidence level, and the dark grey contour represents the 2a confidence level 
around the best-fit. The dotted line denotes matter density Q,Q m (l + z) 3 , and 
the dashed horizontal line denotes ACDM. 



Table 2. The weighted average w (eq 10) over specified redshift ranges for analysis 
using the 'Gold' sample of SNe from [27]. The best-fit value and la deviations from 
the best-fit are shown. 

W 

Q 0m Az = 0-0.414 Az = 0.414-1 Az = 1 - 1.755 
~0^2 -0.847±8-8i§ -O.HSlgl!? O0891H1 
0.3 -1.05318Z -0.1591S 0.1181^1? 

o.4 -1.31011??° -0.2101838 0.21518:815 



for higher values of Oom, the dark energy equation of state has a more negative value 
at present and shows a sharper evolution over redshift. It is interesting to note that 
a constant equation of state which is < — 1 during the evolution of the universe from 
z ~ 1.7 to z = appears not to be supported by the recent SNe data. 

To summarize, our results clearly demonstrate that evolving DE is by no means 
excluded by the most recent SNe observations. On the contrary, our results for pde and 
u>de obtained using the 'Gold' sample of [27] (figures 5, 6) are very similar to the results 
which we obtain using the SNe samples of [21, 22]. For 0.2 < f2 0m < 0.4 the best-fit DE 
model evolves from w < — latz^Otow^Oatz^lin agreement with the results 
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(a) (b) (c) 

Figure 6. The evolution of w(z) with redshift for different values of f2om, 
for the 'Gold' sample of SNe from [27]. The reconstruction is done using the 
polynomial fit to dark energy, equation (5). In each panel, the thick solid line 
shows the best-fit, the light grey contour represents the la confidence level, and 
the dark grey contour represents the 2a confidence level around the best-fit. 
The dashed line represents ACDM. 



of Paper I. 

2.2. DE reconstruction using SNe('Gold')+CMB : 

Observations of the cosmic microwave background and Type la supernovae provide us 
with complementary insight into the nature of dark energy [36, 37, 38, 39]. We may use 
the WMAP result of R = It dz/h(z) = 1.710 ± 0.137 (from WMAP data alone) 

in conjunction with the SNe 'Gold' sample to reconstruct DE. For this purpose, we use 
Q b h 2 = 0.024 and Q 0m h 2 = 0.14 ± 0.02 [3]. To calculate z is we use a fitting function 
given in [40]: 

z ls = 1048[1 + 0.00124(fi 6 /iT a738 ][l + gii^omh 2 ) 92 } , (11) 

where 

9l = 0.078(lM 2 r a238 [l + 39.5(lW) a76 T\ (12) 
( ?2 = 0.56[l + 21.1(fi 5 /i 2 ) 1 - 81 ]- 1 . (13) 

Figure 7 shows our reconstruction of DE obtained with the ansatz (5) using the 
WMAP result together with the SNe 'Gold' sample. The best fit values for this 
reconstruction are: fiom = 0.385, A\ = —2.87,^2 = 1.01, h = 0.60. The best fit dark 
energy density is Pde(^) = 2.475 — 2.87a; + 1.01a; 2 . Note that the best fit p DE decreases 
monotonically with redshift up to z ~ 0.4 and then begins to increase. This reflects the 
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Figure 7. Results from analysis of SNe(Gold)+CMB data, using ansatz (5). 
Q,0m and h are fixed at best-fit values of Oom = 0.385, h = 0.60. Panel (a) 
shows the confidence levels in the (A\,A2) parameter space. The star marks 
the best-fit and the filled circle marks the ACDM point. Panel (b) shows the 
logarithmic variation of dark energy density with redshift. Panel (c) shows the 
evolution of dark energy equation of state with redshift. Panel (d) shows the 
variation of the deceleration parameter with redshift. In all three panels (b), 
(c) and (d), the thick solid line represents the best-fit, the light grey contours 
represent the la confidence level, and the dark grey contours represent the 2a 
confidence levels. The dashed line in panels (b), (c) and (d) represents ACDM, 
and the dotted line in panel (b) represents the matter density. The horizontal 
thick solid line in (d) represents the la limits on the transition redshift at which 
the universe starts accelerating. 
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Table 3. The weighted average w (eq 10) over specified rcdshift ranges for analysis 
from SNe+CMB data. The best-fit value and la deviations from the best-fit are shown. 

W 

fi 0m Az = 0- 0.414 Az = 0.414 - 1 Az = 1 - 1.755 
~0385 =T287™ -0.229±ffi 0.142™ 



phantom nature of DE (w < —1) at lower redshifts. The equation of state behaves as 
before, evolving from wo < — ltoiu~0at;z~l. From the figure 7(d), we see that 
the deceleration parameter q has a value of g = —0.84 ±0.11 at present. The transition 
from deceleration to acceleration (q(z T ) = 0)occurs at a redshift of z T = 0.39 ± 0.03. 
Therefore, from a joint analysis of CMB and SNe data, one may obtain a fairly good 
idea of when the universe began to accelerate. These results demonstrate that the best 
fit to SNe+CMB observations favours evolving DE with a somewhat higher value of f2 0m 
and a slightly lower value of h. Note that the value of f2om is larger than the WMAP 
result f2 0m = 0.27 ± 0.04 obtained using the ACDM prior, but is in agreement with the 
results obtained in [41, 42] (flom — 0.35 ± 0.12) using high redshift clusters. Also note 
that, because of the dust-like behaviour of DE at higher redshifts, the above value of 
fl 0m obtained using equations (2-5) maybe somewhat larger than the values obtained for 
fi m from clustering measurements. The value of h ~ 0.60 obtained for the best fit is in 
tension with the ACDM based result from WMAP, h = 0.73 ±0.03, but is in agreement 
with the observations of [43, 44, 45, 46, 47, 48, 49] which can accommodate lower values 
of h ~ 0.6 (see also the discussion in [50] in this context). We therefore conclude that 
a joint analysis of SNe and CMB data favours evolving DE over ACDM if no priors are 
placed on fiom and h separately (f2o m /i 2 = 0.14 ± 0.02 is assumed). 

It is however a useful exercise to see how the behaviour of DE would change if 
strong priors were imposed on f2 0m and h. We therefore show results obtained by using 
the ACDM based CMB priors [3]: Q 0m = 0.27±0.04 and h = 0.71 ±0.06. The best fit in 
this case has f2 0m = 0.29 and a dark energy density of p DE = 1.23 — 0.81a; ± 0.29a; 2 . The 
equation of state at present is w = —1.10, and it slowly evolves to w(z = 1.75) ~ —0.4. 
The deceleration parameter has a value of go — —0.63 ±0.12 at present and the redshift 
at which the universe begins to accelerate is zt = 0.57 ± 0.07. Figure 8 demonstrates 
that the time evolution of DE is extremely weak in this case and is in good agreement 
with ACDM cosmology (see also [51]). 

3. Conclusions 

In conclusion, we find that the case for evolving dark energy (originally demonstrated in 
Paper I) is upheld by the new supernova data if no priors are imposed on f2 0m and h. For 
a reasonable range of 0.2 < Oom < 0.4, the equation of state of dark energy evolves from 
w < — 1 today to w ~ at z ~ 1. The above result remains in place if we add CMB 
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Figure 8. Results for analysis of SNe(Gold)+CMB data with ACDM-based 
WMAP priors of n 0m = 0.27 ± 0.04 and h = 0.71 ± 0.06, using ansatz (5). 
Panel (a) shows the confidence levels in the (A\,A2) space. The star marks 
the best-fit and the filled circle marks the ACDM point. Panel (b) shows the 
logarithmic variation of dark energy density with redshift. Panel (c) shows the 
evolution of dark energy equation of state with redshift. Panel (d) shows the 
variation of the deceleration parameter with redshift. In all three panels (b), 
(c) and (d), the thick solid line represents the best-fit, the light grey contours 
represent the la confidence level, and the dark grey contours represent the 2a 
confidence levels. The dashed line in panels (b), (c) and (d) represents ACDM, 
and the dotted line in panel (b) represents the matter density. The horizontal 
thick solid line in (d) represents the la limits on the transition redshift at which 
the universe starts accelerating. 
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priors to the analysis. In this case, evolving dark energy with f2 0m — 0.385 and h ~ 0.6 
is favoured over ACDM, and the epoch at which the universe began to accelerate is 
zt = 0.39 ±0.03 within la. However, if we assume strong priors on Oo m and h using the 
ACDM based WMAP results, then the best-fit chooses an Q 0m = 0.29. The evolution 
in the equation of state becomes weaker and is in much better agreement with ACDM. 
The redshift of transition from deceleration to acceleration is z T = 0.57 ±0.07, which is 
closer to the ACDM value of z T ~ 0.7 for this value of f2 0m . Finally, one must note that 
the DE evolution becomes weaker or stronger depending on the subsampling of the SNe 
dataset. A larger number of supernovae at high redshifts, as well as better knowledge 
of the values of H and Q 0m are therefore required before firm conclusions are drawn 
about the nature of dark energy. 
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